function [x_next, y_next, z_next] = Logistic(x, y, z, x_prev, z_prev, r, b)
%QUANTUMLOGISTICMAP 量子 Logistic 系统单步迭代
%   [x_next, y_next, z_next] = quantumLogisticMap(x, y, z, x_prev, z_prev, r, b)
%   输入：
%     x, y, z       — 当前三维状态
%     x_prev, z_prev— 上一时刻 x 和 z（用于差分近似）
%     r             — 分岔参数 γ
%     b             — 参数 β
%   输出：
%     x_next, y_next, z_next — 下一步状态

    % 差分近似
    xdot = x - x_prev;
    zdot = z - z_prev;

    % 映射方程 (Eq.1)
    x_next = r*(x - x^2) - r*y;
    y_next = -y*exp(-2*b) + exp(-b)*r*((2 - x - xdot)*y - x*zdot - xdot*z);
    z_next = -z*exp(-2*b) + exp(-b)*r*(2*(1 - xdot*z)*z - 2*x*y - x);
end
